Derealization in coupled one-dimensional chains 
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A weakly disordered quasi-one-dimensional tight-binding hopping model with N rows is consid- 
ered. The probability distribution of the Landauer conductance is calculated exactly in the middle 
of the band, e = 0, and it is shown that a derealization transition at this energy takes place if and 
only if N is odd. This even-odd effect is explained by level repulsion of the transmission eigenvalues. 

PACS numbers: 72.15.Rn, 11.30. R 



The existence of derealization transitions in a disor- 
dered one-dimensional system is surprising, as it goes 
against the general wisdom that disordered systems in 
less than two dimensions are localized Nevertheless 
a derealization transition in one dimension goes back 
to Dyson's work on models for a glass in 1953 [§,1). 
Dyson's one-dimensional glass is related to a large variety 
of disordered systems: a one-dimensional tight-binding 
model with nearest-neighbor random hopping || , a two- 
dimensional asymmetric random bond Ising model 
which is equivalent to the one-dimensional random quan- 
tum Ising chain || , one-dimensional random bond quan- 
tum XY models || and more generally random XYZ 
spin- 1/2 Heisenbcrg models [Q, and narrow gap semi- 
conductors ||. These models are of current interest in 
view of their rich physics: new universality classes, loga- 
rithmic scaling, the existence of strong fluctuations call- 
ing for a distinction between average and typical prop- 
erties. They might also be useful laboratories to address 
the problem of disorder induced quantum phase transi- 
tions in higher dimensions such as the plateau transition 
between insulating Hall states in the quantum Hall effect 

IS 

The one-dimensional nearest-neighbor random hop- 
ping model is described by the Hamiltonian 



^ — y ] {^n c n c n+l d~ t n C n+l c n 



(1) 



where the operators and c n are creation and annihila- 
tion operators for spinless fermions, respectively, and the 
hopping parameter t n = t + St n consists of a non-random 
part t and a fluctuating part St n . The fundamental sym- 
metry of the Hamiltonian (Q) that distinguishes it from 
one-dimensional systems with on-site disorder is the pres- 
ence of a sublattice (or chiral) symmetry: particles can 
hop only from even to odd- numbered sites. The energy 
e = is special since it corresponds to a logarithmically 
diverging mean density of states |^|. Furthermore, there 
are several independent correlation lengths that diverge 
for e — > indicating that the energy e — represents 
a (disorder induced) quantum critical point |^,^|,0J^]. In- 
particular, at e = the conductance exhibits large fluctu- 



ations superimposed on an algebraically decaying mean 
value Jic[ |. By contrast, for nonzero energy the system 
described by Eq. (Q) is non-criticial resulting in standard 
localized behaviour: A typical sample is well character- 
ized by (logg), which is proportional to L and has rela- 
tively small sample- to-sample fluctuations. 

A different type of derealization in one-dimensional 
disordered systems was considered recently by Hatano 
and Nelson jll| , who considered a chain with on-site dis- 
order and an imaginary vector potential. As a function of 
the strength of the imaginary vector potential, the system 
reaches a critical point and goes through a derealization 
transition. 

The discussion so far applies to the case of strictly one- 
dimensional systems. In this Letter we address the ques- 
tion of whether aspects of the behaviour described above 
carry over to the multi-channel case. Surprisingly, it will 
turn out that the answer depends on the parity of the 
channel number N: For N even, the system behaves very 
much like standard disordered multi-channel wires, i.e. in 
the limit L — > oo all states are localized. However, for N 
odd, precisely one mode remains critical and, moreover, 
exhibits much of the behaviour of the single critical mode 
of strictly one-dimensional systems. For large L, where 
the contribution of all other — localized — modes is ne- 
glegible, the phenomenology of the wire is determined 
by the contribution of the single critical mode, and, in 
this sense, remains critical. To our knowledge, this par- 
ity effect was first noticed by Miller and Wang in their 
study of random flux and passive advection models [ p"2"| . 
However, in that work, the effect has been washed out by 
taking a two-dimensional thermodynamic limit. Keeping 
N finite we here focus on a different regime, where parity 
has pronounced phcnomcnological consequences. 

To probe the onset of critical behaviour as e — > we 
calculate the probability distribution of the conductance. 
As will be shown below the even/odd effect manifests it- 
self in the level repulsion between transmission eigenval- 
ues. We also discuss the effect of staggering in the non- 
random part of the hopping parameter t (connected e.g. 
to a Peierls instability) and establish a relation between 
derealization transition in random hopping models and 
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in non-Hermitian quantum mechanics. 

To be specific, we consider the Hamiltonian 



N 



t n ,ij C n,j C n+l,i 4~ ^n,ij C n+l,i C n,j 



(2) 



where the indices i and j label the N chains. Weak stag- 
gering in the hopping is introduced by setting t n ij = 
tSij + (—l) n t'Sij + 5t n ,ij, where t' -C i. We distinguish 
between the cases in which time reversal symmetry is 
present {(3 = 1, t n ,ij real) from those where it is absent 
(ft = 2, t n .ij complex). The weakly fluctuating parts of 
the hopping amplitudes 8t ni ij are taken to be indepen- 
dent and Gaussian distributed, with zero mean and with 
variance (St n ij St* l ij ) = ftv 2 /j, where 7 = ftN + 2-/3. 

Upon linearization of the spectrum in the vicinity 
of the Fermi energy e = 0, the lattice model can 
be approximated by a continuum model obeying the 
Schrodinger equation 

N 

ei>i{y) = ^2hij(y)4>j{y), i=l,...,N, (3a) 

hij = iv F S i: jcridy + Vij(y)ai + Wij(y)a 2 . (3b) 

Here ip is a two-component wavefunction, correspond- 
ing to even and odd-numbered sites in the original lat- 
tice model, y — 2na, a being the lattice constant, and 
vf = 2ta is the Fermi velocity. The sublattice symme- 
try of the lattice model (^|) translates to 03/1^ 03 = —hij, 
which we refer to as chiral symmetry. The chiral symme- 
try distinguishes this system from one-dimensional sys- 
tems with on-site disorder, which do not show a dereal- 
ization transition. The random potentials v and w are 
Hermitian (vij = v* t , Wij = w*i) while, in the presence 
of time-reversal symmetry, one has the further condition 
= —v*j and Wij = u>,* ■. Apart from the symme- 
try constraints, the random potentials are independent 
and Gaussian distributed, with mean (vij(y)) = and 
(wij(y)) = 2t'5ij, and variance (v 2 = 2v 2 a(3^f~ 1 ) 

(Sv ij (y)Sv ij (y')*) = v 2 S(y - y'){l - S P i5 lj ) 
(Swij(y)Swij(y')*) = v 2 S(y - y')(l + 5p\6ij). 

In order to find the conductance at zero energy, we 
calculate the distribution of the 2N x 2N transfer matrix 
M , which relates wavefunctions at the left and right of 
a disordered strip of length L Q. The eigenvalues of 
MM\ which arise in inverse pairs exp(±2xj), determine 
the transmission eigenvalues Tj = 1/ cosh 2 Xj and hence 
the conductance g through the Landauer formula 
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(4) 



In the absence of disorder, all exponents Xj are zero, and 
conduction is perfect, g = N. On the other hand, trans- 
mission is exponentially suppressed if all larger 



than unity. The Xj S clX6 related to the channel-dependent 
localization lengths = L/\xj\. The largest length £ de- 
termines the exponential decay of the conductance g and 
serves as the localization length of the total system of 
coupled chains. 

To compute the distribution of M, we use the Fokker- 
Planck approach pioneered for disordered wires with 
random on-site disorder by Dorokhov G3| and Mello, 
Pereyra, and Kumar fl5|| . Following the method of Refs. 
|l4],[l5| , we first consider the case of disorder confined to 
a small strip < y < SL. Denoting the wave function for 
y < by ipj(L) and for y > SL by ipj(R), we find 



N 



(5) 



fc=i 



where the (random) transfer matrix M of the slice reads 



SL 



M = T y exp I Vp 1 I dy [iv(y) - a 3 w{y)} 



(6) 



Here T y denotes the ordering operator for the y- 
integration. 

For any given realization of the disorder, the transfer 
matrix has the following symmetry properties: 



(73 Mas — M (chiral symmetry), 
Ma\M^ = o\ (flux conservation), 
M* = M (time reversal). 



(7a) 
(7b) 
(7c) 



Taking the symmetries (]?]) into account, we find that the 
transfer matrix can be parameterized as 



M = u exp(xa3)v , 



(8) 



where u and v are the tensor product of N x N unitary 
matrices (orthogonal if ft = 1) with the 2x2 unit ma- 
trix and x is a diagonal N x N matrix with real diagonal 
elements x\, . . . , xn- The numbers x±, . . . , xn are the 
radial coordinates of the transfer matrix (eigenvalues of 
ilogMAfT), and the matrices u and v are the angular 
coordinates. In contrast to systems without chiral sym- 
metry, the Xj can be both positive and negative. 

The transfer matrix of a system of length L is found 
by multiplication of the transfer matrices of the many 
individual slices of width SL. As each multiplication re- 
sults in a small change of the radial coordinates xj , they 
perform a "Brownian motion" ]l4| , |T5| |. Upon multiplica- 
tion with the transfer matrix of a slice of width SL, we 
find that the radial coordinates Xj change according to 
Xj — > Xj + Sxj, where the first two moments of the in- 
crement Sxj , averaged over the disorder configuration in 
the added slice, are given by 
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(*■':, >al = 2^ | ' ' X! coth(ij •'•/, : ) • (9a) 



(5xj5x k ) SL = j^tjk- 



(9b) 



Here the mean free path £ and dimensionless staggering- 
disorder ratio / read 



I = v jr /4v A a, f = jt'v F /v 2 a/3 



(9c) 



The first term on the r.h.s. of Eq. ( |9a| ) results in a si- 
multaneous drift of all radial coordinates Xj. The sec- 
ond term describes repulsion between nearby Xj in the 
Brownian motion process. The Fokker-Planck equation 
corresponding to Eq. (||) reads 



„dP 



N 



d 



J = Y\ I smh(xj 

k>j 



Xk)\ 



dxj 



(10a) 
(10b) 



The initial condition corresponding to perfect transmis- 
sion at L = is P(xi, . . . ,xn; 0) = Y[j 3( x j)- 

The Fokker-Planck equation (|o|) is the central result 
of this Letter |lq| . It contains all information on the 
transport properties of the random hopping system at 
zero energy. Eq. ([h]) is the chiral analogue of the so- 
called Dorokhov-Mello-Pereyra-Kumar (DMPK) equa- 
tion, which governs the evolution of the transmission 
eigenvalues of a disordered wire @-||. The key differ- 
ence between the two equations is the presence of "mirror 
imaged" eigenvalues Xj in the DMPK equation, which are 
absent in Eq. (fTo[). [For wires with on-site disorder, the 
eigenvalues xjnot only repel from different eigenvalues 
Xk, c.f. Eq. (pa), but also from the "mirror image" —Xk', 
in particular, Xj and — Xj repel.] 

In the absence of time- reversal symmetry ((3 = 2), the 
DMPK equation has been solved exactly by Beenakker 
and Rejaei pj^ by a mapping to a problem of non- 
interacting fcrmions. Using the method of Ref. (l7|, we 
have been able to find an exact solution of Eq. ( |To[ ) for 
(3 = 2. It reads 



P = C (L)JJexp [-fx. 



,N£ 

T 



X JJ (x k - xj) smh(x k - xj), 
k>j 



(11) 



where c(L) is a normalization constant. The exact so- 
lution (|l0|) has a formal analogy to the distribution of 
eigenvalues of a random matrix: it consists of a pair inter- 
action and a potential part. However, while for random 
matrices the eigenvalue interaction is quadratic, here we 
find a more complicated level repulsion. Comparing our 
result ( |TT| ) to the exact solution of Beenakker and Rejaei, 



we note the absence of the mirror-image eigenvalues in 
the interaction and potential factors. No exact solution 
of Eq. (|h]) for (3 = 1 could be found. 

To determine whether the system is at a critical point, 
we investigate the distribution of Xj's for L — > oo, which 
can be obtained from Eq. (^ for both (3 — 1 and (3 = 2. 
For L 3> N£, the radial coordinates Xj are well sepa- 
rated, say x\ <C ... <^ xn- We then find from Eq. (||) 
that the "dynamics" of the Xj's (j = 1, . . . , N) separate, 
and that they show small Gaussian fluctuations around 
equidistant equilibrium positions, 

(xj) =(N+l-2j- f)L(3/2l 1 , vaxxj = L/jL (12) 

This is the so-called "crystallization of transmission 
eigenvalues" which is a signature of localization in 
wires with on-site disorder. Transmission is exponen- 
tially suppressed if all radial coordinates larger 
than unity, c.f. Eq. (^). For on-site disorder all Xj's grow 
linearly with L |L3| , which inevitably leads to strong lo- 
calization. [Within the framework of the DMPK equa- 
tion, this results from the repulsion between Xj and the 
mirror image — Xj.] The situation is different for the cou- 
pled random hopping chains, where we find from Eq. ( |T^ ) 
that the radial coordinate Xj remains (on average) close 
to zero, thus resulting in a delocalized state and a critical 
point, provided 



N + 1 - 2j - / = 0. 



(13) 



As a result, in the absence of staggering (/ = 0), a crit- 
ical point exists only if the number of chains is odd. If 
there is no staggering, an even number of coupled ran- 
dom hopping chains show an exponential decay of the 
conductance. The conductance distribution at the criti- 
cal point follows directly from the Landauer formula (Q) 
and the Gaussian distribution of the radial coordinate 
Xj. As fluctuations of Xj around zero are large [see Eq. 
(O)], the conductance at the critical point shows large 
sample-to-sample fluctuations, and the random hopping 
chains at the critical point can by no means be regarded 
as a "good conductor". 

The parity effect for the presence of a critical point 
in the absence of staggering can be understood from the 
"level repulsion" of the variables Xj. In the large-L limit, 
where x\ <C • • ■ <C xn, the coordinates Xj repel by con- 
stant forces, see Eq. ([)a]). For an even number of chan- 
nels, there is a net force on all Xj's, driving them away 
from and resulting in an exponential suppression of the 
conductance (see Fig. |l[a). However, as is depicted in Fig. 
[l]b, if the number of channels is odd, there is no force on 
the middle exponent a;(jv+i)/2- Therefore, this variable 
will remain close to zero and give rise to a diverging lo- 
calization length £ = L/|a;(jv+i)/2l an d a critical state. 
For comparison, in the case of a wire with on-site disor- 
der, the repulsion between Xj and its mirror image — Xj 



3 



results in a nonvanishing force for all radial coordinates 
@ (see Fig. |c). 

By fine tuning the staggering parameter / ( pcj ) , which 
measures the ratio of the uniform staggering if and the 
random disorder strength v, an additional [N/2] critical 
points can be reached, both for even and odd number of 
chains ( N/2] is the largest integer < N/2). According 
to Eq. (1J), as the staggering parameter / approaches 
the critical value / = 2j — N — 1, the localization length 
£ = £j = L/\xj\ diverges with (critical) exponent 1. 
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FIG. 1. The parity effect results from level repulsion be- 
tween the transmission eigenvalues, (a) For an even number 
of chains, all radial coordinates Xj are repelled away from 0, 
while (b) for an odd number of chains a;(jv+i)/2 remains close 
to 0. (c) Repulsion from mirror images for a wire with on-site 
disorder results in a positive driving force for all Xj. 

The fact that we find exponential suppression of the 
conductance if Eq. ( |i~3| ) is not obeyed, may be due to ei- 
ther the existence of localized states, or to a gap in the 
spectrum. For instance, in Peierls materials, staggering 
opens a gap in the excitation spectrum, which explains 
the exponential suppression of the transmission even for 
zero disorder. The presence of disorder leads to a finite 
density of states below the excitation gap, but these sub- 
gap states are localized except for critical realizations of 
the disorder strength that satisfy Eq. (13). 

To close, we discuss the relation between the critical 
points for the multi-chain random hopping model stud- 
ied in this Letter and the Hatano-Nelson delocalization 
transition infone-dimensional) non-Hermitian quantum 
mechanics . The relation between the two systems is 
established through the "method of Hermitization" jl8| , 
in which the non-Hermitian problem with "Hamiltonian" 
h at (complex) energy z is made Hermitian by consider- 
ing the Hamiltonian H z = uiRe (h — z) + o^Im (h — z). 
An eigenfunction of h at eigenvalue z is an eigenfunc- 
tion of H z at eigenvalue and vice versa. For complex 
disorder in the the non-Hermitian system, we find that 
the iV-chain non-Hcrmitian problem maps to 2N coupled 
chains with Hermitian quantum mechanics and with chi- 
ral symmetry. The staggering parameter in the chiral 
system plays the role of an imaginary vector potential 
considered by Hatano and Nelson |11|. Thus, compar- 
ing the non-hermitian problem with the random hopping 
chain containing an even number of rows, we deduce 



that, in the absence of the imaginary vector potential, 
the non-Hermitian system is localized. Since the imag- 
inary vector potential maps to the staggering, a series 
of N critical points (and the corresponding branches of 
delocalized states with complex energy, see Ref . [jllj ) can 
then be obtained by tuning the values of the imaginary 
vector potential. 
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